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Abstract 



Dense neutron matter with recently developed BSkl9 and BSk21 Skyrme effective forces is considered in magnetic fields up to 
£S) 10 20 G at zero temperature. The breaking of the rotational symmetry by the magnetic field leads to the differentiation between the 
q pressures along and perpendicular to the field direction which becomes significant in the fields H > H t h ~ 10 18 G. The longitudinal 
<D pressure vanishes in the critical field 10 18 < H c < 10 19 G, resulting in the longitudinal instability of neutron matter. For the Skyrme 
£^ force fitted to the stiffer underlying equation of state (BSk21 vs. BSkl9) the threshold H,h and critical H c magnetic fields become 
1 larger. The longitudinal and transverse pressures as well as the anisotropic equation of state of neutron matter are determined under 
the conditions relevant for the cores of magnetars. 

. . Keywords: Neutron matter, strong magnetic field, magnetar model, spin polarization, pressure anisotropy, Fermi liquid approach. 
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1. Introduction 

Magnetars are strongly magnetized neutron stars 01 with 
emissions powered by the dissipation of magnetic energy. Ac- 
cording to one of the conjectures, magnetars can be the source 
of the extremely powerful short-duration y-ray bursts J2, S 0. 
13] . The magnetic field strength at the surface of a magnetar is 
of about 10 14 -10 15 G 00]. Such huge magnetic fields can be 
inferred from observations of magnetar periods and spin-down 
rates, or from hydrogen spectral lines. In the interior of a mag- 
netar the magnetic field strength may be even larger, reaching 
values of about 10 18 G H>@]. Under such circumstances, the is- 
sue of interest is the behavior of neutron star matter in a strong 
magnetic field fill HE} El El El El. 
■ In the recent study 111 111 , neutron star matter was approx- 
imated by pure neutron matter in a model with the effec- 
tive nuclear forces. It was shown that the behavior of spin 
polarization of neutron matter in the high density region in 
a strong magnetic field crucially depends on whether neu- 
tron matter develops spontaneous spin polarization (in the ab- 
sence of a magnetic field) at several times nuclear matter sat- 
uration density, or the appearance of spontaneous polariza- 
tion is not allowed at the relevant densities (or delayed to 
much higher densities). The first case is usual for the Skyrme 
forces OS El [H [H M, El EHH Hill IH, while the sec- 
ond one is characteristic for the realistic nucleon-nucleon (NN) 
interaction 1 27 , 28 , H 3^, 31 , 32 , 33 , 34 ] . In the former case, a 
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ferromagnetic transition to a totally spin polarized state occurs 
while in the latter case a ferromagnetic transition is excluded 
at all relevant densities and spin polarization remains quite low 
even in the high density region. 

The scenario for the evolution of spin polarization at high 
densities in which the spontaneous ferromagnetic transition 
in neutron matter is absent was considered for the magnetic 
fields up to 10 18 G lllll . However, it was argued in the re- 
cent study 13511 that in the core of a magnetar the local values 
of the magnetic field strength could be as large as 10 20 G, if 
to assume the inhomogeneous distribution of the matter density 
and magnetic field inside a neutron star, or to allow the forma- 
tion of a quark core in the high-density interior of a neutron 
star (concerning the last point, see also Ref. l36ll ). Under such 
circumstances, a different scenario is possible in which a field- 
induced ferromagnetic phase transition occurs in the magnetar 
core. This idea was explored in the recent research l37ll . where 
it was shown within the framework of a lowest constrained vari- 
ational approach with the Argonne NN potential that a fully 
spin polarized state in neutron matter could be formed in the 
magnetic field H > 10 19 G. Note, however, that, as was pointed 
out in the works IB5LI38I1 . in such ultrastrong magnetic fields the 
breaking of the 0(3) rotational symmetry by the magnetic field 
results in the anisotropy of the total pressure, having a smaller 
value along than perpendicular to the field direction. The pos- 
sible outcome could be the gravitational collapse of a magnetar 
along the magnetic field, if the magnetic field strength is large 
enough. Thus, exploring the possibility of a field-induced fer- 
romagnetic phase transition in neutron matter in a strong mag- 
netic field, the effect of the pressure anisotropy has to be taken 
into account because this kind of instability could prevent the 
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formation of a fully polarized state in neutron matter. This ef- 
fect was not considered in Ref. 13711 . thus, leaving the possibil- 
ity of the formation of a fully polarized state of neutron spins 
in a strong magnetic field open. In the given study, we pro- 
vide a fully self-consistent calculation of the thermodynamic 
quantities of spin polarized neutron matter taking into account 
the appearance of the pressure anisotropy in a strong magnetic 
field. We consider spin polarization phenomena in a degenerate 
magnetized system of strongly interacting neutrons within the 

HHHm, unlike 



framework of a Fermi liquid approach ||3S 
to the previous works 1351 13811 . where interparticle interactions 
were switched off. 

Note that recently new parametrizations of Skyrme forces 
were suggested, BSkl9-BSk21 M43I1 . aimed to avoid the spon- 
taneous spin instability of nuclear matter at densities beyond 
the nuclear saturation density for vanishing temperature. This 
is achieved by adding new density-dependent terms to the stan- 
dard Skyrme interaction. The BSkl9 parametrization was con- 
strained to reproduce the equation of state (EoS) of nonpolar- 
ized neutron matter [44] obtained in variational calculation with 
the use of the realistic Urbana vu nucleon-nucleon potential 
and the three-body force called there TNI. The BSk20 force cor- 
responds to the stiffer EoS l45ll . obtained in variational calcula- 
tion with the use of the realistic Argonne Vi% two-body potential 
and the semiphenomenological UIX* three-body force which 
includes also a relativistic boost correction 5v. Even a stiffer 
neutron matter EoS was suggested in the Brueckner-Hartree- 
Fock calculation of Ref. 14611 based on the same two-body 
potential and a more realistic three-body force containing dif- 
ferent meson-exchange contributions. This EoS is the under- 
lying one for the BSk21 Skyrme interaction. Further we would 
like to contrast the results obtained with the Skyrme forces con- 



strained to soft li44ll and most stiff [;46|] underlying EoS, and, by 
this reason, choose the BSkl9 and BSk21 parametrizations in 
the subsequent analysis. 

At this point, it is worthy to note that we consider thermody- 
namic properties of spin polarized states in neutron matter in a 
strong magnetic field up to the high density region relevant for 
astrophysics. Nevertheless, we take into account the nucleon 
degrees of freedom only, although other degrees of freedom, 
such as pions, hyperons, kaons, or quarks could be important at 
such high densities. 

2. Basic equations 

The normal (nonsuperfluid) states of neutron matter are de- 
scribed by the normal distribution function of neutrons f KlKl - 
Tr ga^a Kl , where k = (p, cr), p is momentum, cr is the projec- 
tion of spin on the third axis, and g is the density matrix of the 
system 124 l25l 12611 . The energy of the system is specified as a 
functional of the distribution function f,E- E(f), and deter- 
mines the single particle energy 

dE(f) 
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df K2 



(1) 



The self-consistent matrix equation for determining the distri- 
bution function / follows from the minimum condition of the 



thermodynamic potential [39, 4jJ and is 

/ = {exp(T e + T, ■ fi„cn + Y 4 ) + 1} 
= {exp(F £ + l}" 1 • 



(2) 



Here the quantities e, T, and Y4 are matrices in the space of k 
variables, with (Y iA ) KlK2 = Y iA 6 KlK2 , Y = l/T, Yi = -Hi /T and 
Y4 = -fiu/T being the Lagrange multipliers, jUo being the chem- 
ical potential of neutrons, and T the temperature. In Eq. (0, 
H„ = -1.9130427(5W ~ -6.031 ■ 10~ 18 MeV/G is the neutron 
magnetic moment 14711 (fi^ being the nuclear magneton), cr, are 
the Pauli matrices. Note that, unlike to Refs. 014111511 . the term 
with the external magnetic field H is not included in the single 
particle energy s but is separately introduced in the exponent of 
the Fermi distribution ©. 

Further it will be assumed that the third axis is directed along 
the external magnetic field H. Given the possibility for align- 
ment of neutron spins along or opposite to the magnetic field 
H, the normal distribution function of neutrons and the matrix 
quantity £ (which we will also call a single particle energy) can 
be expanded in the Pauli matrices <x, in spin space 



/(P) = /o(P)o"0 + / 3 (p)cr 3 , 
£(P) =£o(p)cr +£3(p)cr 3 . 



(3) 
(4) 



The distribution functions /o,/" 3 satisfy the normalization 
conditions 
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7y Zjf^)^ei-ei = Ap. 



(5) 
(6) 



Here g — g^ + qi is the total density of neutron matter, £>f and gi 
are the neutron number densities with spin up and spin down, 
respectively. The quantity Ag may be regarded as the neutron 
spin order parameter which determines the magnetization of the 
system M = fi„Ag. The spin ordering of neutrons can also be 
characterized by the spin polarization parameter 



n 



Ag 
Q ' 



The magnetization may contribute to the internal magnetic field 
B = H + 4-nM. However, we will assume, analogously to the 
previous studies ll9l [Tll[l4ll . that, because of the tiny value of the 
neutron magnetic moment, the contribution of the magnetiza- 
tion to the inner magnetic field B remains small for all relevant 
densities and magnetic field strengths, and, hence, 



BsH. 



(7) 



Indeed, e.g., the magnetic field necessary to produce a fully 
polarized spin state is, at least, greater than 10 19 G at the densi- 
ties relevant for the cores of magnetars (as will be shown later), 
while for totally spin polarized neutron matter with the density 
g — 1 frrT 3 the contribution of the term with the magnetization 
to the inner magnetic field amounts at AnM ^ 1.2 • 10 17 G. 
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In order to get the self-consistent equations for the compo- 
nents of the single particle energy, one has to set the energy 
functional of the system. It represents the sum of the matter 
and field energy contributions 



E{f,H) = E m (f) + ^- f V. 

tin 



(8) 



The matter energy is the sum of the kinetic and Fermi-liquid 
interaction energy terms II251 12611 



EJJ) = E Q (f) + E int (f), 

£o(/)=2^£ (p)/o(p), 
p 

Eintif) = 2&>(P)/0(P) + £3(P)/ 3 (P)}, 



(9) 



where 



q 



(10) 

(ii) 



Here e (p) = ^- is the free single particle spectrum, mo is 
the bare mass of a neutron, t/ '(k), i/"(k) are the normal Fermi 
liquid (FL) amplitudes, and eg, £3 are the FL corrections to the 
free single particle spectrum. Taking into account Eqs. ([T]),© 
and (0, expressions for the components of the single particle 
energy read 

ft(p) = e (p) + eo(p) - fi , 6(p) = + (12) 

In Eqs. (fT2l i. the quantities eg, £3 are the functionals of the 
distribution functions /c/3 which, using Eqs. (f2| and (0, can 
be expressed, in turn, through the quantities 



fa = ^{n(^+) + n(a)-)}, 
fi = ^{«(w+) - «(w_)}, 



(13) 
(14) 



where 



n((o ± ) = {exp(F w ± ) + 1) 



The quantity <x> ± , being the exponent in the Fermi distribu- 
tion function n, plays the role of the quasiparticle spectrum. 
The branches a> ± correspond to neutrons with spin up and spin 
down. 

Thus, Eqs. (fT2l-([T4l> form the self-consistency equations for 
the components of the single particle energy, which should be 
solved jointly with the normalization conditions (0, ©. 

The pressures (longitudinal and transverse with respect to the 
direction of the magnetic field) in the system are related to the 
diagonal elements of the stress tensor whose explicit expression 
reads |48 | 



CTik = 



Sit + 



An 



Here 



?-f H 



(15) 



(16) 



f H = 7?(E - TS) - HM is the Helmholtz free energy density. 
For the isotropic medium, the stress tensor (JT3J is symmetric. 

The transverse p, and longitudinal pi pressures are deter- 
mined from the formulas 

Pt = -0"11 = -0"22. Pi = -0"33- 

At zero temperature, using Eqs. ((H), $15[ , one can get the ap- 
proximate expressions 



dg 

/de m \ H L 



8tt' 

2 



(17) 
(18) 



where e m is the matter energy density, and we disregarded the 
higher order small terms containing M. The structure of the 
pressures p t and pi is different that reflects the breaking of 
the rotational symmetry by the magnetic field. In ultrastrong 
magnetic fields, the quadratic on the magnetic field term (the 
Maxwell term) will be dominating, leading to increasing the 
transverse pressure and to decreasing the longitudinal pressure. 
Hence, at some critical magnetic field, the longitudinal pressure 
will vanish, resulting in the longitudinal instability of neutron 
matter. The question then arises: What is the magnitude of the 
critical field and the corresponding maximum degree of spin 
polarization in neutron matter? 

3. Longitudinal and transverse pressures. Anisotropic EoS 
at zero temperature 

In order to solve the self-consistent equations, we utilize 
the BSkl9 and BSk21 parametrizations of the Skyrme interac- 
tion, developed in Ref. 14311 and generalizing the conventional 
Skyrme parametrizations. By choosing these Skyrme forces, 
we would like to study the influence of the underlying EoS, 
to which these Skyrme forces were constrained, on thermody- 
namic quantities of strongly magnetized dense neutron matter. 

The normal FL amplitudes in Eqs. dTOb.dTTb can be related to 
the parameters of the Skyrme interaction by formulas [49] 

I/£(k) = 2f (l - x ) + 1^(1 - x 3 ) + -|[f!(l - xi) (19) 



+ f 4 (l - x^g 8 + 3/ 2 (l + x 2 ) + 3/ 5 (l + x 5 )£> r ]k\ 
1 = -2/ (l - x ) - |e"(l - x 3 ) + |j[f2(l + x 2 ) 

+ f 5 (l + X 5 )Q? -h(l- JCi) - f 4 (l - X 4 )(f]k 2 . 



In these equations, the terms with the factors t\ and fg are the 
additional density-dependent terms generalizing the t\ and t 2 
terms in the conventional form of the Skyrme interaction f50ll . 
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Figure 1 : (Color online) Neutron spin polarization parameter as a function of 
the magnetic field H for the Skyrme forces BSkl9 (solid lines) and BSk21 
(dashed lines) at zero temperature and fixed values of the density, g = 4po (two 
lower branches) and g = 6go (two upper branches). The vertical arrows indicate 
the maximum magnitude of spin polarization attainable for the corresponding 
Skyrme force at the given density, see further details in the text. 



They were added to the usual form with the aim to avoid the 
appearance of spontaneous spin instabilities in nuclear and neu- 
tron matter at high densities. Specific values of the parameters 
tj, Xj, a, ft and y as well as of the nuclear saturation density go 
for each parametrization are given in Ref. 114311 . 

Now we present the results of the numerical solution of the 
self-consistency equations. Fig. 1 shows the spin polarization 
parameter of neutron matter as a function of the magnetic field 
H at two different values of the neutron matter density, g = 4go 
and g = 6go, which can be relevant for the central regions of 
a magnetar. It is seen that the impact of the magnetic field 
remains small up to the field strength 10 17 G. For the BSk21 
force (stiff underlying EoS), the magnitude of the spin polar- 
ization parameter is smaller that that for the BSkl9 force (soft 
underlying EoS). For both parametrizations, the larger the den- 
sity is, the smaller the effect produced by the magnetic field 
on spin polarization of neutron matter. At the magnetic field 
H = 10 18 G, usually considered as the maximum magnetic field 
strength in the core of a magnetar (according to a scalar virial 
theorem 115 lh . the magnitude of the spin polarization param- 
eter doesn't exceed 25% for the BSkl9 force and 8% for the 
BSk21 force (for the densities under consideration). However, 
the situation changes if the larger magnetic fields are allowable: 
With further increasing the magnetic field strength, the magni- 
tude of the spin polarization parameter increases till it reaches 
the limiting value IT = — 1, corresponding to a fully spin polar- 
ized state. For example, this happens at H a 1.3 • 10 19 G for 
the BSkl9 force and at H * 4.3 ■ 10 19 G for the BSk21 force at 
g = 4po, i-e-, certainly, for magnetic fields H > 10 19 G. Never- 
theless, we should check whether the formation of a fully spin 
polarized state in a strong magnetic field is actually possible 
by calculating the anisotropic pressure in dense neutron matter. 
The meaning of the vertical arrows in Fig. 1 is explained later 
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Figure 2: (Color online) (a) Pressures, longitudinal (descending branches) and 
transverse (ascending branches), as functions of the magnetic field H for the 
Skyrme forces BSkl9 (solid lines) and BSk21 (dashed lines) at zero tempera- 
ture and fixed values of the density, g = Ago (two lower branching curves) and 
g = 6go (two upper branching curves), (b) Same as in the top panel but for the 
normalized difference between the transverse and longitudinal pressures. The 
vertical arrows in the lower panel indicate the points corresponding to the onset 
of the longitudinal instability in neutron matter. 



in the text. 

Fig. 2a shows the pressures (longitudinal and transverse) 
in neutron matter as functions of the magnetic field H at the 
same densities, g — Ago and g — 6go. The upper branches 
in the branching curves correspond to the transverse pressure, 
the lower ones to the longitudinal pressure. First, it is clearly 
seen that up to some threshold magnetic field the difference be- 
tween transverse and longitudinal pressures is unessential that 
corresponds to the isotropic regime. Beyond this threshold 
magnetic field strength, the anisotropic regime holds for which 
the transverse pressure increases with H while the longitudinal 
pressure decreases. The stiffer the underlying EoS is (BSk21 
vs. BSkl9), the larger the pressure, transverse p, or longitu- 
dinal pi. Also, the increase of the density has the same ef- 
fect on the pressures p, and pi as stiffening of the underlying 
EoS. The most important feature is that the longitudinal pres- 
sure vanishes at some critical magnetic field H c marking the 
onset of the longitudinal collapse of a neutron star. For exam- 
ple, H c ~ 1.7 10 18 G for BSkl9 force and // f * 3.2 10 18 Gfor 
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BSk21 force at g = 4g , and H c « 3.4 ■ 10 18 G for BSkl9 force 
and « 6.3 • 10 18 G for BSk21 force at g = 6g . In all cases 
under consideration, this critical value doesn't exceed 10 19 G. 



In Ref. 13511 . the critical field for a relativistic dense gas of free 
charged fermions was found to be close to 10 19 G. 

The magnitude of the spin polarization parameter II cannot 
also exceed some limiting value corresponding to the critical 
field H c . These maximum values of the ITs magnitude are 
shown in Fig. 1 by the vertical arrows. In particular, II f as -0.42 
for BSkl9 force and II, * -0.26 for BSk21 force at g = 4g , 
and II C as -0.13 for BSkl9 force and II, as -0.20 for BSk21 
force at g = 6go. As can be inferred from these values, the ap- 
pearance of the negative longitudinal pressure in an ultrastrong 
magnetic field prevents the formation of a fully spin polarized 
state in the core of a magnetar. Therefore, only the onset of a 
field-induced ferromagnetic phase transition, or its close vicin- 
ity, can be catched under increasing the magnetic field strength 
in dense neutron matter. A complete spin polarization in the 
magnetar core is not allowed by the appearance of the negative 
pressure along the direction of the magnetic field, contrary to 
the conclusion of Ref. 13711 where the pressure anisotropy in a 
strong magnetic field was disregarded. 

Fig. 2b shows the difference between the transverse and lon- 
gitudinal pressures normalized to the value of the pressure po in 
the isotropic regime (which corresponds to the weak field limit 
with pi = p, = p ): 

o Pt - Pi 



Po 

It is quite reasonable to admit that, when the anisotropic regime 
sets in, the splitting between the transverse and longitudinal 
pressures becomes comparable with the value of the pressure 
in the isotropic regime 13 511 . Applying for the transition from 
the isotropic regime to the anisotropic one the approximate 
criterion 6-1, the transition occurs at the threshold field 
H th ~ 1.2 • 10 18 G for BSkl9 force and H, h as 2.2 • 10 18 G 
for BSk21 force at g = 4g , and at H th as 2.3 ■ 10 18 G for BSkl9 
force and H th as 4.6 • 10 18 G for BSk21 force at g = 6g . In all 
cases under consideration, the threshold field H,i, is larger than 
10 18 G, and, hence, the isotropic regime holds for the fields up 
to 10 18 G. For comparison, the threshold field for a relativis- 
tic dense gas of free charged fermions was found to be about 
10 17 G 0511 (without including the anomalous magnetic mo- 
ments of fermions). For a degenerate gas of free neutrons the 
model dependent estimate gives H,i, 4.5 • 10 18 G 0~8ll (in- 
cluding the neutron anomalous magnetic moment). The nor- 
malized splitting of the transverse and longitudinal pressures 
increases more rapidly with the magnetic field at the smaller 
density and/or for the Skyrme force BSkl9 with the softer un- 
derlying EoS. The vertical arrows in Fig. 2b indicate the points 
corresponding to the onset of the longitudinal instability in neu- 
tron matter. Since the threshold field H,i, is less than the critical 
field H c for the appearance of the longitudinal instability, the 
anisotropic regime can be relevant for the core of a magnetar. 
The maximum allowable normalized splitting of the pressures 
corresponding to the critical field H c is 6 ~ 2. If the anisotropic 
regime sets in, a neutron star has the oblate form. Thus, as fol- 
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Figure 3: (Color online) The energy density of the system as a function of: 
(a) the transverse pressure p,, (b) the longitudinal pressure p; for the Skyrme 
forces BSkl9 (solid lines) and BSk21 (dashed lines) at zero temperature and 
fixed values of the density, g = 4po and g = 6go. The meaning of the vertical 
arrows in the top panel is the same as in Fig. 2. In the bottom panel, the physical 
region corresponds to p; > 0. 



lows from the preceding discussions, in the anisotropic regime 
the pressure anisotropy plays an important role in determining 
the spin structure and configuration of a neutron star. 

Because of the pressure anisotropy, the EoS of neutron mat- 
ter in a strong magnetic field is also anisotropic. Fig. 3 shows 
the dependence of the energy density of the system on the trans- 
verse pressure (top panel) and on the longitudinal pressure (bot- 
tom panel) at the same densities considered above. Since in an 
ultrastrong magnetic field the dominant Maxwell term enters 
the pressure p, and energy density with positive sign and the 
pressure pi with negative sign, the energy density is the increas- 
ing function of p, and decreasing function of pi. In the case of 
e{p t ) dependence, at the given density, the same p, corresponds 
to the larger magnetic field H for the BSkl9 force compared 
with the BSk21 force (see Fig. 2a). The overall effect of two 
factors (the stiffness/softness of the underlying EoS and mag- 
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netic field) will be the larger value of the energy density at the 
given p t and density for the BSkl9 force compared with the 
BSk21 force (see Fig. 3a). The analogous arguments show that, 
for the given Skyrme parametrization and at the given p,, the 
energy density is larger for the smaller density. In the case of 
e{pi) dependence, at the given density, the same pi corresponds 
to the smaller magnetic field H for the BSkl9 force compared 
with the BSk21 force (see Fig. 2a). Hence, the energy density 
at the given pi and density is larger for the BSk21 force than 
that for the BSkl9 force (see Fig. 3b). Analogously, for the 
given Skyrme parametrization and at the given pi, the energy 
density is larger for the larger density. In the bottom panel, the 
physical region corresponds to the positive values of the longi- 
tudinal pressure. Note that, because of the validity of the ap- 
proximation given by Eq. (01, the energy density containing the 
field energy contribution is practically indistinguishable from 
the Helmholtz free energy density. 

It is worthy to notice that the occurrence of the longitudinal 
instability in a strong magnetic field will lead to the compres- 
sion of a neutron star along the magnetic field with a subsequent 
increase of the density. Such an increase can eventually cause 
the appearance of new particle species. Already for the condi- 
tions of high density considered in this study one can assume 
that a deconfined phase of quarks exists in the interior of a neu- 
tron star. Then a hybrid star consisting of deconfined quark 
matter in the core and nuclear matter in the outer layers can be 
regarded as a relevant astrophysical object. As follows from 
the general arguments presented in our study, the critical mag- 
netic field strength at which the longitudinal pressure vanishes 
inevitably exists for a hybrid star as well. The determination of 
the corresponding critical value needs a separate investigation. 
As a possible guess for this value, one can refer to the results 
of the study of quark matter within the MIT bag model l35ll . 
where it was estimated to be close to 10 19 G (with the bag con- 
stant put to zero in the final expressions). Also, the mass-radius 
relationship is a relevant characteristic of neutron stars. Usually 
it is found by solving the Tolman-Oppenheimer-Volkoff (TOV) 
equations |52{] for a spherically symmetric and static neutron 
star. In an ultrastrong magnetic field, the EoS becomes essen- 
tially anisotropic. Unlike to the standard scheme, the mass- 
radius relationship should be found by the self-consistent treat- 
ment of the anisotropic EoS and axisymmetric TOV equations 
substituting the conventional TOV equations in the case of an 
axisymmetric neutron star. 

Note that in this research we have studied the impact of a 
strong magnetic field on thermodynamic properties of dense 
neutron matter at zero temperature. It would be also of inter- 
est to extend this research to finite temperatures relevant for 
proto-neutron stars which can lead to a number of interesting 
effects, such as, e.g., an unusual behavior of the entropy of a 



spin polarized state 153115411 . 

In summary, we have considered spin polarized states in 
dense neutron matter in the model with the Skyrme effective 
NN interaction (BSkl9 and BSk21 parametrizations) under the 
presence of strong magnetic fields up to 10 20 G. It has been 
shown that in the magnetic field H > H,i, ~ 10 18 G the pressure 
anisotropy has a significant impact on thermodynamic proper- 



ties of neutron matter. In particular, vanishing of the pressure 
along the direction of the magnetic field in the critical field 
H c > H t h leads to the appearance of the longitudinal instability 
in neutron matter. For the Skyrme force with the stiffer under- 
lying EoS (BSk21 vs. BSkl9), the threshold H, h and critical 
H c magnetic fields become larger. The increase of the density 
of neutron matter also leads to increasing the fields //,/, and H c . 
Even in the extreme scenario with the most stiff underlying EoS 
considered in this work and at the densities about several times 
nuclear saturation density, the critical field H c doesn't exceed 
10 19 G which can be considered as the upper bound on the mag- 
netic field strength inside a magnetar. Our calculations show 
that the appearance of the longitudinal instability prevents the 
formation of a fully spin polarized state in neutron matter, and 
only the states with mild spin polarization can be developed. 
The longitudinal and transverse pressures and anisotropic EoS 
of neutron matter in a strong magnetic field have been deter- 
mined at the densities relevant for the cores of magnetars. The 
obtained results can be of importance in the structure studies of 
strongly magnetized neutron stars. 
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